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Superfluid phenomena can be explained in terms of the topologies of the order parameter and of 
the confining vessel. For example, currents in a toroidal vessel can be characterized by a discrete and 
conserved quantity, the winding number. In trapped Bose-Einstein condensates, the topology of 
the trap can be characterized by the topology of the Thomas-Fermi surface of its N-particle ground 
state. This can be altered during an experiment, so that a toroidal trap may deform into a more 
spherical shape, allowing an initially persistent current to decay into singly-quantized vortices. We 
investigate such a procedure numerically, and confirm that the Thomas-Fermi prescription for the 
trap topology gives an accurate picture of vortex formation. 



The intriguing beauty of the concept of topological 
quantum numbers has ever since its introduction into 
physics drawn great attention to states that are topo- 
logically nontrivial Q. In a system that has undergone 
a second order phase transition these numbers charac- 
terize the order parameter field. For a Bose-condensed 
system the order parameter is the phase of the macro- 
scopic wave-function, and vortex line defects, as well as 
persistent currents in multiply-connected vessels, are ex- 
amples of topologically nontrivial states that have been 
well investigated in superfluid "^He and ■^He Q. Indeed, 
these classic superfluids are sufficiently strongly interact- 
ing systems that the topology of the order parameter 
field becomes especially important, because it is insensi- 
tive to the microphysical details for which we lack simple 
and reliable theories. Dilute samples of Bose-condensed 
alkali atoms however, are weakly interacting, and 
the numerically tractable Gross-Pitaevskii (GP) mean 
field theory describes them with high quantitative accu- 
racy. They are also amenable to a wide range of exper- 
imental controls, including the time-dependent manipu- 
lation of the inhomogeneous external trapping fields. In 
this paper we propose an experimental approach which 
takes advantage of this capability to probe the continu- 
ous range of physics between vortices and persistent cur- 
rents, by adiabatically deforming a trapping potential so 
that the topology it imposes on the condensate effectively 
changes. This causes an initially uniform current to break 
up, at some point, into a number of vortices. The tech- 
nique thus offers an adiabatic complement to the stirring 
experiments on critical velocity in condensates Q. 

The topological quantum number of a vortex defect or 
a persistent current can be found by writing the wave- 
function of the superfluid as the product of a modulus 
and a phase, 'J = \fpe^^ ■ The condition of superfluidity 
means that the density field p is strongly fixed energet- 
ically, so that the remaining low-energy degrees of free- 
dom are those of the current velocity field, given by the 
gradient of the phase, v — {Ti/m)V9, with m being the 
particle mass. Because 9 must be single-valued modulo 




FIG. 1: Density of (a) a persistent current in a toroidal trap 
with ht = 25 and cr^ = 5 (see eq. (^)) and (b) a vortex state 
in a harmonic trap. For both plots n = 1. 



27r, the circulation calculated by integrating v along an 
arbitrary closed contour C, 
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must be 2'K{fi/m)K. The integer k will be unchanged by 
any continuous deformation of the field 9: it is a topolog- 
ical invariant known as the winding number. If the entire 
superfluid sample has no 'holes', however, then the con- 
tour C can be shrunk to zero size, so that n is changed 
to zero by a continuous process, after all. The paradox is 
resolved by the fact that the winding number can change 
if p vanishes at any point on the contour, because there 
9 becomes undefined. Hence a simply connected sample 
with K ^ must contain such topological defects - vortex 
cores. A multiply connected sample can have a winding 
number without defects, in which case the current cannot 
decay - unless vortices nucleate at the sample surface, 
and then drift across the current. Thus the interplay 
between topological defects in the order parameter, and 
the topology of the confining vessel, gives a qualitative 
explanation of superfluidity, which is a necessary prelim- 
inary to any microscopic considerations (of how vortices 
nucleate, etc.). 

For liquid helium there is a sharp distinction between 
topological defects, which are microscopically small, and 
the topology of the sample's macroscopic container. For 



dilute Bose-Einstein condensates, in contrast, a vortex 
core may be only a couple of orders of magnitude smaller 
than the entire sample. So if we are to use trap topol- 
ogy to get a simple understanding of how cold atomic 
currents can be stable, and how they can decay, then 
we must be precise about how to distinguish the 'intrin- 
sic topology' of a condensate, due to vortices, from the 
topology which is imposed by the trap. Since topology 
characterizes closed surfaces, trap topology must be the 
topology of an equipotential surface V{r) = Vq ■ For 
traps of finite strength, this topology will in general de- 
pend on the choice of Vq. For a given number of trapped 
atoms N, however, the natural candidate is the ground 
state Thomas-Fermi (TF) surface V{r) = fi, where fJ.{N) 
is the condensate chemical potential. It is a main result 
of this paper to show that this prescription does indeed 
provide a very useful measure of the effect of the trap 
on vortex nucleation in circulating condensates. The dis- 
tinction between a condensate with a vortex, and a con- 
densate circulating in a toroidal trap, remains from one 
point of view only a matter of degree (see Fig. |l|); but this 
definition of trap topology captures the fact that there is 
a physically important difference of scale. 

We simulate the evolution of a dilute Bose-Einstein 
condensate held initially in a toroidal trap 0, ob- 

tained by focusing a blue-detuned laser beam into the 
center of a rotationally symmetric harmonic trap ||, ^ . 
For numerical tractability we consider the limit of a 
quasi-two-dimensional trap, and neglect the third di- 
mension entirely. Assuming the profile of the laser to 
be Gaussian, the resulting potential Vt(r) can be written 
as 

Vt{r) = K^ + hte-^ , (2) 

where ht and a determine the height and width of the 
central peak. The subscript t indicates that the height 
is to be modified adiabatically during the evolution [[lO[ . 
Using natural units we scale the length in units of the 
harmonic oscillator ground state size oq — y^h/muj and 
the energy in units of the trap frequency huj. Because 
the potential is changed adiabatically, the state of the 
system at time t is assumed to be determined by the 
instantaneous time-independent GP equation 

/^tVt--^A2i/'t-f FtVt+slV'tlVt , (3) 

where the non-linear coupling constant in the quasi- 
2D trap is given by <? = ina/ao, with a represent- 
ing the atomic s-wave scattering length. We calcu- 
late the sequence of instantaneously stationary states 
of the adiabatically evolving system by propagating the 
time-dependent GP equation in imaginary time , and 
changing the potential slightly in every time step. We be- 
gin with an initial ipo with k > 1, with ho, (Tq such that 
this state is (locally) stable. As we slowly lower ht, we 
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FIG. 2: Breaking up of a k = 2 current, showing the respec- 
tive density (upper row) and phase (lower row) distributions. 
The height ht is lowered from (a) to (c), with a constant, 
yielding (a) the stable k = 2 current, then (6) a locally stable 
state with two k = 1 vortices, and finally (c) the non-rotating 
ground state. 

find in all cases that at some critical height he instabil- 
ity occurs. In Fig. |^ we show a generic example of this 
effect, in which a current with k = 2 breaks up into two 
vortices. (The rotational symmetry is broken through nu- 
merical noise, with vortices forming at angles unrelated 
to the spatial grid.) 

Our goal is to investigate quantitatively the behaviour 
of this critical height he as a function of k and a. Al- 
though the imaginary time solution of eq. (|^) does give 
the correct evolution of the system, the vortex cores form 
at such low densities that our numerical resolution is in- 
sufficient to determine the precise barrier height at which 
they first form. We overcome this limitation by numer- 
ically solving the imaginary time Bogoliubov equations 
around the solutions of the GP equation, using the ansatz 

^it^ipt + ^ , (4) 
so that the imaginary time Bogoliubov equation reads 

(/i + £)(p = - i (^drr + ^dr + ^^00^ Lp 

+ Vir)^ + 2U\^t\''^ + U\i;t\^^*- (5) 

We solve Eq. (||) for the lowest eigenvalue e by relaxation 
in imaginary time, using a higher order Runge-Kutta 
method. The eigenfunctions ip found by this procedure 
have no direct physical significance; but the first appear- 
ance of a negative s indicates precisely when instability 
occurs. The nature of the instability is then determined 
by the imaginary time GP evolution as outlined above. 
We repeat these calculations for a range of k and a; the 
results are shown in Fig. |[ 

The most salient feature of the results is a 'floor' effect: 
in many cases, with wider barriers and higher winding 
numbers, the critical barrier height is equal to the chem- 
ical potential. The critical barrier height may rise above 
this 'floor', if the barrier is narrow enough or the winding 
number small enough. These two basic results suggest 
one natural interpretation: the spawning of the vortices 



FIG. 3: Stability and instability as a function of k and . 
States 'inside' the bars are unstable, he increases with k and 
decreases with a The chemical potential increases from « 18 
to ^ ~ 27 with increasing k and is independent on the scale 
of the plot of a. Fields where ht = Q (large k, small a) were 
beyond our numerical capabilities because the steepness of 
the central peak necessary to obtain a stationary states could 
not be resolved sufficiently with reasonable spatial grid size. 
Here g\ip\^ = 1000. 



is governed by two factors, one of which is the topology 
of the ground state TF surface. We show below that 
these two factors are conceptually worth distinguishing, 
but are not actually physically independent. 

The first factor in vortex formation is instability of sur- 
face excitations [|l2|. If the velocity of the atoms at the 
inner border exceeds the local critical velocity, Vc, vor- 
tices form while the sample is still multiply connected on 
the trap scale. Lowering the central barrier shrinks the 
inner radius of the toroidal TF surface, and with fixed 
circulation this raises the velocity at this inner radius. 
If K is high enough, stability may fail when the inner 
toroidal radius is so large that the curvature of the inner 
TF surface is irrelevant, and we are effectively seeing the 
critical velocity for flow past a straight TF surface. In 
this limit we can expect to recover the critical velocity 
for flow perpendicular to a linear ramp potential (because 
the vortices will nucleate within a narrow boundary layer, 
within which the potential is effectively linear). From di- 
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mensional analysis it then follows that Vc oc {TiF/w?^ 
where F is the derivative of the trapping potential, at 
the TF surface, in the direction normal to that surface 
[ p2[ . A comparison of this approximation with the re- 
sults from the numerical solution of eq. (^ is shown in 
Fig. ^. With appropriately chosen proportionality con- 
stants, which evidently reflect the role of curvature, the 
agreement is excellent. Since the scaling of the critical 
velocity with the surface force persists quite well in al- 
most all of our cases, it is evident that the picture of 
vortices forming through surface instability is at least 
qualitatively valid even in cases where the concept of a 



FIG. 4: Critical velocity of the current at the inner TF surface 
as a function of the width of the central potential. Circles: 
numerical solutions of eq. (|^; Full lines: TF approximation 
as described in the text with the proportionality constants 
A{n) = [2, 2.5, 2.9, 3.3, 3.8, 3.9]. Note the different scale in the 
plot for Hi = 1. 



surface may well be suspect (such as those in which the 
condensate density is very low where the vortices form). 

We can supplement our understanding of vortex forma- 
tion in this doubtful regime, and explain the floor effect 
in Fig. ^ by invoking trap topology as our second fac- 
tor. Lowering ht below the chemical potential abruptly 
changes the trap topology, and turns a persistent cur- 
rent of winding number |k| > 1 into a multiply quantized 
vortex. As is well known, a multiply quantized vortex is 
unstable against breaking up into unit vortices, because 
the energy of a single vortex is (x k^. Of course, this topo- 
logical factor is not a separate mechanism, in addition to 
the microphysical mechanisms of instability at supercrit- 
ical velocity. It should be obvious that as the inner TF 
surface disappears, the velocity at that surface becomes 
arbitrarily large. But trap topology is a concise sum- 
mary of microphysical details, which can quickly identify 
a whole regime in which unit vortices will appear. 

On the other hand, there can be qualitative features of 
vortex behaviour which depend on more specific details 
of the trapping potential than its topology. Our results 
corroborate previous observations that instability to vor- 
tex formation need not mean that the resulting vortices 
penetrate the current iQ. Fig. ||(b) is not an interme- 
diate stage in a real-time process, but a local energy 
minimum for a particular potential height. It shows that 
vortices may be stably trapped near the inner TF surface. 
From the point of view of superfluidity, this indicates 
that a persistent current may be nonlinearly stabilized 
despite being linearly unstable. (That is, the linear in- 
stability can saturate nonlinearly without destruction of 
the current.) From the point of view of vortex dynamics, 
this trapping feature would allow study of interactions 
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between vortices. Since the energy of a vortex scales 
with the density of the background fluid, vortices tend 
to 'float' towards the potential barrier iQ; but the well- 
known Coulomb-like repulsion between vortices circulat- 
ing in the same direction tends, in the circular geometry, 
to drive them radially outwards. A too gentle slope of the 
central potential, or too many vortices, allows this repul- 
sion to dominate buoyancy, and the vortices are expelled 
through the outer TF surface. As a result the winding 
number of the current circulating in the trap is decreased 
by the number of exiting vortices. So one corollary appli- 
cation of our topology-changing procedure is stabilizing 
vortices against drifting out of the sample. 

The tunable central barrier can also be useful for creat- 
ing circulating currents in the first place. Simulations in- 
dicate that phase imprinting fl^ methods should be par- 
ticularly enhanced, because the problematic core region 
of the mask becomes irrelevant. Phase imprinting should 
also be effective in even less trivial topologies: with two 
Gaussian barriers, and a slightly more complicated mask, 
one can make vortex pairs, or vortex-anti-vortex pairs. 
And once quantized circulation has been established 
around the barriers, they may be moved independently. 
Finally, the non-equilibrium phenomenon of spontaneous 
creation of winding number in rapid quenches has been 
predicted to be more readily observable in toroidal traps 
[ p^t , and with a tunable barrier one could assess this de- 
pendence quantitatively. 

In summary, therefore, we have proposed and exam- 
ined a method of effectively changing the topology of a 
trap confining a condensate with a persistent current. We 
have computed critical parameters for vortex formation, 
and identified the effect of trap topology. The proposed 
technique would allow the precisely controlled produc- 
tion of vortices, and may help to shed new light on the 
general role of topology in superfluid phenomena. 
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